Temperature dependence of polarization relaxation in semiconductor quantum dots 
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The decay time of the linear polarization degree of the luminescence in strongly confined semicon- 
ductor quantum dots with asymmetrical shape is calculated in the frame of second-order quasielastic 
interaction between quantum dot charge carriers and LO phonons. The phonon bottleneck does not 
prevent significantly the relaxation processes and the calculated decay times can be of the order of 
a few tens picoseconds at temperature T ~ 100K, consistent with recent experiments by Paillard et 
al. [Phys. Rev. Lett. 86, 1634 (2001)]. 



PACS numbers: 78.67.Hc, 72.25.Rb, 63.20.Ls 

The discrete nature of the energy spectrum in semi- 
conductor quantum dots (QDs) is supposed to lead to 
a strong suppression of spin relaxation 1 where promis- 
ing applications for new spin dependent electronic de- 
vices have been predicted 2 ' 3 . The current interest in ma- 
nipulating semiconductor spins for applications is based 
on the ability to control and maintain spin coherence 
over practical length and time scales. Optical pump- 
ing experiments have indeed given good indications of 
a slowing down of the carrier spin relaxation processes 
in QD compared to bulk or quantum well structures. 4,5 
Recently, a detailed time-resolved investigation of the in- 
trinsic spin dynamics in InAs/GaAs QDs under strictly 
resonant excitation has been reported by Paillard et al. 6 
They demonstrated that at low temperature the carrier 
spins are totally frozen on the scale of the exciton life- 
time. A rapid temporal decay of the linear polarization 
degree is, however, found above 30K. 

As a possible intrinsic mechanism for temporal decay 
of the linear polarization, we propose here the second- 
order quasielastic interaction between QD carriers and 
LO phonons as sketched in Fig. 1(b). A similar mecha- 
nism was suggested by Uskov et al. 7 to explain the tem- 
perature dependence of a broadening of the zero-phonon 
line. We will show that such relaxation processes are not 
suppressed (strongly) by the "phonon bottleneck" and 
can lead to a decay time of the linear polarization of 
some tens of picoseconds at ~ 100K. 

Let us first recall that self-organized QDs are usually 
strained and have an asymmetrical shape with a height 
smaller than the base size. The upper valence band 
in such QDs with the zinc-blende lattice is split into a 
heavy-hole band with the angular momentum projections 
jh z = ±3/2 and a light-hole band with jh,z = ±1/2 at 
the center of the Brillouin-zone (here the growth direc- 
tion z is chosen as the quantization axis). The conduc- 
tion band is s-like, with two spin states s e ^ z = ±1/2. The 
heavy-hole exciton quartet is characterized by the total 
angular momentum projections J z — ±1, ±2. The ra- 
diative states J z = ±1 and nonradiative ones, J z = ±2, 
are split by the electron-hole exchange interaction (so- 



called singlet-triplet splitting A s t). 8-10 In the case of an 
asymmetric confinement potential in the plane of QDs 
the symmetry of the system is lowered and the exchange 
interaction is no longer isotropic. 1,11 As a result, both 
doublets are split into singlets, as is shown schemati- 
cally in Fig. 1. The radiative doublet | ± 1) is split by 
an anisotropic exchange into the states labeled \X) = 
(|l)+t|-l))/V2and|Y) = (|1)-»|-1))/V2, giving rise to 
optical transitions which are linearly polarized along the 
[110] and [110] directions, respectively. 12,13 Continuous 
wave single dot spectroscopy experiments have clearly ev- 
idenced these two linearly polarized lines in self-organized 
InGaAs QDs. 13 An analogous splitting of the exciton 
states has been observed in studies of single QD's formed 
at GaAs/AlGaAs interfaces. 14 The observation of the 
"optical orientation-optical alignment" or "alignment- 
orientation" transformation in a magnetic field makes 
it possible to determine the magnitude of the splitting 
without resolving the fine-structure spectrally. 6,12,15 The 
anisotropic exchange splitting originates from the elon- 
gation of the QDs. 1 ' 12,16 The calculated and measured 
magnitudes of this splitting reach some tens or even hun- 
dreds of /LteV. 6, 12-14 
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FIG. 1. Schematic diagram showing (a) the sublevels of the 
exciton ground state and (b) the second-order phonon scatter- 
ing process within the radiative doublet via the first excited 
state. Optically-inactive states are shown by dashed lines. 



1 



In the case of (pulsed) resonant excitation which is 
linearly polarized, say, along the Y-direction (<j y - exci- 
tation), an exciton in the Y-sublevel is created. If there 
is no relaxation to the X-sublevel, the following emission 
occurs from the same Y- sublevel, as is shown in Fig. 1(a), 
and the linear polarization degree of the luminescence 
Plin = (I Y -I X )/(I X +I Y ) (where I x and I Y denote the 
X and Y linearly polarized luminescence components) re- 
mains constant within the exciton lifetime, Pu n = 1. The 
possible relaxation process is shown schematically in Fig. 
1 (b) . This is the second-order exciton-phonon scattering 
process in which one phonon is absorbed and another 
one is emitted. The scattering events occur via the ex- 
cited states of the exciton (with negligible exchange split- 
ting) which can couple to both X- and F-exciton states. 
These processes lead to the exciton transitions between 
the ground Y— and X- sublevels and result in the tem- 
poral dependence of Pu n . The (complete) rate equations 
governing the populations fy and fx of the Y- and X- 
sublevels read 
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where t sc is the scattering time for an exciton relax- 
ing from the YpQ-state to the X(Y)-state and T ra d 
is the radiative recombination time. For the relevant 
initial conditions fy(0) = fo and /x(0) = 0, Eqs. (1) 
lead to an exponential decay of the polarization de- 
gree Pun = (f Y - fx)/{fv + fx) = exp(-t/T po i) with 

As was noted above, in experiments by Paillard et 
al. 6 no decay of Pu n is observed on the exciton lifetime 
scale at T < 30K, but the decay time of Pu n drops from 
~ 3.5ns at 40K down to ~ 50 ps at 80K. In order to ex- 
plain these findings we consider the contribution of LO 
phonons to the second-order scattering processes above. 
In this case a strong (nonlinear) temperature dependence 
of a relaxation time t sc can be expected in the tempera- 
ture range ksT smaller than the LO phonon energy HQo- 

Our calculation of the polarization decay time is based 
on the Fermi golden rule for second order processes. To 
simplify matters, we restrict ourselves to a certain range 
of parameters where the following conditions hold: 
(i) We assume that the phonon energy Mlo is in the vicin- 
ity of the excitation energy E from the exciton ground 
state doublet (labelled "0" ) to the first excited state (la- 
belled "1") so that contributions of higher excited states 
can be neglected in the first approximation. As a result, 
we have 

1 _ 2tt ^ I Mf Mf Mf Mf 2 

q,k 

xN ? (N^ + l) <5(A exc -7% + MV), (2) 
where k, fig, T and = l/{e m z/ kBT -1) arc the phonon 



wave vector, frequency, linewidth, and thermal distribu- 
tion function, respectively 

(ii) To evaluate the exciton-phonon matrix elements 
M2 1 , we restrict ourselves to the Frohlich interaction 17 
of bulk LO phonons with carriers in strongly confined 
flat QDs (strained in the growth direction). For a strong 
confinement, the exciton states are reasonably well ap- 
proximated by noninteracting electron-hole pair states. 
In addition, the electron mass is much smaller than the 
heavy-hole mass and the first excited exciton state is 
mainly composed of a ground state electron and an ex- 
cited hole. As a result, the electron does not couple to 
phonons because of the orthogonality of the initial and 
final hole states. Thus, M2 1 reduces to the matrix ele- 
ment constructed with the hole envelope wave functions 

and *f(r*), 

= ( (s^W/ 2 { ^ (f)|e , * W(f)>j (3) 

qVVK 

where q = \q\, V is the normalization volume, and 
1/k = 1/eoo — l/ £ o ( £ o an d £00 denote the static and 
high-frequency limits of the dielectric function). 

(iii) Equation (3) is mostly affected by the oscillatory 
behavior of the exponential function rather than by the 
detailed shape of the wave functions. Therefore, we may 
furthermore neglect a weak lateral anisotropy of the QD 
and, for simplicity, consider a parabolic potential of the 
form ~ ajj(x 2 + y 2 ) + a^z 2 . The potential strengths a\\ 
and a z are inversely proportional to the base length Ly 
and a height L z of the QD, respectively. In this case the 
hole envelope wave functions are given by 

*oVl = ^ exp{-5V +y 2 ) fz 2 }, 

= \/2a,| x <(f). (4) 

In addition, we have neglected the heavy-hole-light-hole 
mixing because in the case of strained flat QDs the en- 
ergy separation between the first heavy-hole and light 
hole levels is large compared to the lateral confinement 
energies. 

(iv) It follows from Eqs. (3,4) that only LO-phonons 
with small wave numbers q < ot\\ <~ L^ 1 contribute 
substantially to the scattering rate Eq. (2). Therefore, 
the phonon dispersion can be ignored in the Bose func- 
tions, N^,N^ — > 7V . In the delta-function, however, 
the dispersion is essential and it will be approximated 
by an isotropic relation, hfl^ = TiOo — vk 2 ■ Below, we 
restrict our consideration to energy window 0.57iQo > 
\E - M7 | > 0.05^o ~ 1.5meV (hfl ~ 30meV). In 
this case \Eq — %Clo\ is larger than the exchange splitting 
A exc , the effective width of the phonon band va 2 , as well 
as the phonon line width T. 18 Therefore, these contribu- 
tions can be neglected in the energy denominators of Eq. 
(2). 
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FIG. 2. Phonon (a) and level detuning (b) functions. 

After some algebra, we obtain our main result 
1 P(w) 
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F(z) = ern(x/i) (1 + l/2z) - (9) 

where To is the characteristic decay time of the QD 
which depends on the elect ron-phonon coupling strength 
k and the phonon dispersion parameter v. The parame- 
ter w — (Eo — Tiflo)/hrto defines the "phonon detuning" 
between the interlevel spacing in the QD and the LO 
phonon energy, whereas A = A exc /(2va 2 ) defines the 
"mismatch" between the radiative level separation and 
the effective width of the phonon band 19 ' 20 . The param- 
eter = \ —{L z jL\\) 2 measures the anisotropy of the QD 
(for flat QDs under consideration > and close to one) 
and erfi(z) = crf(zz)/z is the imaginary error function 21 . 

From the structure of Eq. (5) we can draw qualitative 
conclusions about the influence of temperature, level and 
phonon dctunings, and size and shape of the QD on the 
polarization decay time. 
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FIG. 3. Polarization decay time as a function of the phonon 
detuning parameter for A = 5 and (3 = 0.75. 



The temperature dependence is determined by the 
phonon number No and r po / decreases exponentially with 
temperature if fc^T <§C KQo. 

The level detuning function for flat QDs with the 
anisotropy parameter — 0.91(L 2 /L|| = 0.3), = 
0.84(L Z /L|| = 0.4), and = 0.75(L Z /L\\ = 0.5) is shown 

in Fig. 2(a). D(A,0) decreases monotonically with A 
and has an exponential limit for A ^> 1. Note that 
D(A, 0) enters in Eq. (5) due to the condition of the 
energy conservation during the scattering process (cfm. 
5- function in Eq. (2)). Furthermore, the level splitting 
A exc plays the role of the energy difference (KQo — Eo) 
when compared with an inelastic first-order scattering 
process. 22 As it is well known, such processes are strongly 
suppressed by the phonon bottleneck effect. The second- 
order processes under consideration are, however, more 
efficient. The point is that an anisotropic exchange split- 
ting A exc <~ O.lmcV is by itself relatively small 6 ' 13 , while 
some specific conditions are required to obtain similar 
magnitudes for the energy difference (HClo — Eo). In ad- 
dition, the parameter A = (A exc /va 2 ) is not expected 
to change very strongly with the lateral size of the QD 
due to the size dependence of the exchange splitting 
A exc . 12,23 Thus, A is determined mostly by the lateral 
anisotropy of the QDs and can be of the order of unity 
(using for A exc and va 2 the values A exc <~ O.lmeV and 

va 2 ~ O.OlmeV, noted above, one obtains that A ~ 5). 
It can be seen in Fig. 2(a) also that the level detun- 
ing function D(A,0) is modified by changing the ratio 
L z /L\\. For instance, it increases by a factor of <~ 1.3 at 
~ 0.91 compared to <~ 0.75 in the case of A = 2, by 
a factor of ~ 1.9 at A = 5, and by a factor of ~ 2.6 at 
A = 7. 
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The size dependence of polarization relaxation is de- 
termined by the phonon detuning function P(w) which 
is shown in Fig. 2(b). The nonexponential dependence 
of P(w) on w results in a relatively weak dependence of 
the polarization decay on the phonon detuning. In ad- 
dition, even at relatively large detunings, e.g., w ~ 0.5 
or w ~ —0.3, P(w) is of the order of unity. The relax- 
ation occurs for both positive (Eo > HClo) an d negative 
(Eq < Hflo) detunings but there is some asymmetry. 

Below we give some numerical estimates of the polar- 
ization decay time in InAs QDs which can be considered 
as typical for III-V compounds. Material parameters 
used are: Ml ~ 30.2meV, e ~ 15.15, ~ 12.25, 24 
and v ~ 0.1(nm) 2 meV 19 which implies To — 0.07ps. For 
T <C 350K the number of phonons No is very small so 
that Tpoi decreases exponentially with temperature, e.g., 
by a factor of ~ 80 between T = 40K and T = 80K which 
is in agreement with experiment. 6 Because the numerical 
values of phonon and level detuning, as well as the QD 
anisotropy, are usually not known exactly from experi- 
ments, w, A and (3 are used as adjustable parameters. As 
an illustration, the decay time of the polarization degree 
is shown in Fig. 3 as a function of the phonon detun- 
ing parameter w at different temperatures. The decay 
time Tp i drops from a few hundred (units) of nanosec- 
onds at T = 30K down to a few tens (and even units) 
of picoseconds at T = 120K, depending on the phonon 
detuning. For w = —0.15 (i.e. E Q - KQo — — 4.5meV), 
A = 7 and (3 = 0.75 which seem to be plausible for 
the experiments by Paillard et al., 6 the calculated decay 
times Tp i ~ 3.5ns (40K) and r po / ~ 44ps (80K) closely 
agree with the experimental results. 

In conclusion, the second-order quasielastic interaction 
between charge carriers and LO phonons in strongly con- 
fined asymmetric QDs is identified as an intrinsic mech- 
anism of the temporal decay of the linear polarization 
degree of the luminescence. Despite of the apparent ob- 
struction by the phonon bottleneck and the level detun- 
ing effects, the proposed mechanism leads to decay times 
of the order of a few tens of picoseconds (or even smaller) 
at temperature T ~ 100K which are in agreement with 
experiments. The relaxation processes are more efficient 
in flat QDs with a weak lateral anisotropy. Equation (5) 
may be useful for the optimation of QD-structures for 
application in spintronic devices. 
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